The P-wave [cs}[cs] tetraquark state: y(4260) or y(4660)? 
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The mass of P-wave cs-scalar-diquark cs-scalar-antidiquark state is computed in the framework 
of QCD sum rules. The result 4.69 ± 0.36 GeV is in good agreement with the experimental value of 
Y(4660) but higher than Y(4260)'s, which supports the P-wave [cs][cs] configuration for Y(4660) 
while disfavors the interpretation of V(4260) as the P-wave [cs][cs] state. In the same picture, the 
mass of P-wave [bs][6s] is predicted to be 11.19 ± 0.49 GeV. 
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PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 



I. INTRODUCTION 

Fruitful heavy hadrons have been observed by far, some of which attribute to the J PC = 1~~ family, e.g. 
Y(4260), Y (4360), and Y (4660). The observation of Y (4260) was first announced by BABAR Collaboration 
[lj, which was confirmed later by both CLEO Collaboration Q and Belle Collaboration A fit to the 
resonance yields a mass 42631* MeV [J. Subsequently, Y(4360) [|J0 and Y(4660) were reported by 
BaBar Collaboration and Belle Collaboration, masses of which are 4361±9±9 MeV and 4664±11±5 MeV, 
respectively. Since then, these states have inspired intensive theoretical speculations. Concre tely , Y(4260) 
is proposed as a hybrid charmonium a molecular state 0, a conventional M>(45 ) llOj . an wxd 
molecular state [11], a A C A C baryonium state [12j, and a D\D or DqD* hadronic molecule [13|j ; Y(4360) is 
interpreted as the candidate of the charmonium hybrid or a 3 3 Z?i cc state 14[; Y_(4660) is suggested to be 



onium state 
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16|, a / (980)*' bound state [17lll8j, a 6 3 Si state [19(, 



a 5 3 S*i charmonium a bary 

and a 5 3 Si — 4 3 P>i mixing state 20(. Besides, many other renewed works |21| have appeared continually. 

In the tetraquark picture, Y(4260) is deciphered as the P-wave [cs][cs] state [22I (named as Yr cs i here), 
however, some authors do not go along with the assumption and figure that Y(4260) cannot be a P- 
wave charm-strange diquark-antidiquark 0. Otherwise, some researchers study Y(4660) as a charm- 
strange tetraquark state 24j . Under such a circumstance, it is interesting and necessary to make clear 
whether Y(4260) can be interpreted as the P-wave [cs][cs] state or F(4660) can be a candidate of the 
Y[ cs ] . Indubitably, the quantitative investigation of Yj cs ] 's mass is very instructive for comprehending its 
structure, but it is quite difficult to extract hadronic spectrum information from the QCD basic theory. 
Fortunately, one can make use of QCD sum rules [3] (for reviews see 3-2^] and references therein), which 
is entrenched in the QCD first principle. Just in this work, we devote to reckon the mass of Y[ cs i through 
the QCD sum rule, to study whether Y(4260) or F(4660) can be a P-wave [cs][cs] state. In addition, 
Yb (10890) [30j, |3l| has been interpreted as a P-wave [6g][6g] tetraquark state 32j. Similarly, the bottom 
counterpart [bs] [bs] for Yj c;5 ] could exist, thereby Yj& s ]'s mass is also predicted here. 

The paper is planned as follows. The QCD sum rule for the tetraquark state is introduced in Sec. [TTJ 
and both the phenomenological and QCD side are derived, followed by the numerical analysis and some 
discussions in Sec. IIIII Section IIVI is a brief summary. 



II. THE P-WAVE [Qs][Qs] QCD SUM RULE 
The QCD sum rule bridges the gap between the hadron phenomenology and the quark-gluon interactions. 



By analogy with the structure of P-wave [Qq][Qq] in Ref. 33 1, the Yiq s i is a J PC = 1~~ bound diquark- 
antidiquark state having the flavor content Y^q^ = [Qs] [Qs] with the spin and orbital momentum numbers: 
S[q s ] = 0, StQgi = 0, SiQ s MQg] = 0, and £[q s i[qs] = 1. For the interpolating current, a derivative could be 
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included in order to generate Lsq^Qs] — 1- Presently, one constructs the tetraquark state current from 
diquark-antidiquark configuration of fields, while constructs the molecular state current from meson-meson 
type of fields. While these two types of currents can be related to each other by Fiertz rearrangements, 
the relations are suppressed by a typical color and Dirac factor so that one could obtain a reliable sum 
rule only if one has chosen the appropriate current to have a maximum overlap with the physical state 



(on this point, there are some calculations and discussions in the XII. Appendix in Ref. [34j ) . Concretely, 
it will have a maximum overlap for the tetraqurk state using the diquark-antidiquark current and the 
sum rule can reproduce the physical mass well, whereas the overlap for the tetraqurk state employing a 
meson-meson type of current will be small and the sum rule will not be able to reproduce the mass well. 
Thus, the following form of current could be constructed for Y^q s ^ 

f = e abc e dec {s T a C l5 Q b )D^- Sdlz CQ T e ). (I) 

Here the index T means matrix transposition, C is the charge conjugation matrix, denotes the covariant 
derivative, as well as a, b, c, d, and e are color indices. 

To derive the mass sum rule, one starts from the two-point correlator 

IT"V) - * / d i xe i ^{0\T\j^x)f+(0)}\0). (2) 

Lorentz covariance implies that the two-point correlator can be generally parameterized as 

ifV) = (^f -snn (1 V) + n«V). (3) 

q A q A 

The part of the correlator proportional to is chosen to attain the sum rule here. Phenomenologically, 
IlW(g 2 ) can be expressed as 



nW( 9 2 ) = J^-L + L / ds lnUi 1 ^ + subtractions, (4) 



[AW] 2 1 f°° Imn( 1 )P hc n (.s) 

where Mh denotes the mass of the hadronic resonance. In the OPE side, IlW(<? 2 ) can be written as 

/■do „OPE M 

nW( 9 2 )=/ ds p - w (5) 

where the spectral density is given by p^^'^'(s) — —Jmll^ (s). After equating the two sides, assuming 
quark-hadron duality, and making a Borel transform, the sum rule can be written as 

[A«] 2 e- M «/ M2 = f° dsp 0PE (s)e- s / M2 . (6) 

J (2m Q +2m s ) 2 

Eliminating the hadronic coupling constant one could yield 

Ml = I" dsp 0PE se-^ 2 / f S ° dsp ° P \-^ 2 . (7) 

J (2m Q +2m B ) 2 J (2m Q +2m B ) 2 

For the OPE calculations, one works at leadingorder in a s and considers condensates up to dimension 
six, with the similar techniques developed in [35|, |36( . The s quark is dealt as a light one and the diagrams 
arc considered up to the order m s . To keep the heavy-quark mass finite, one uses the momentum-space 



expression for the heavy-quark propagator, and the expressions with two and three gluons attached [37 1 
are used. The light-quark part of the correlation function is calculated in the coordinate space and then 
Fourier-transformed to the momentum space in D dimension. The resulting light-quark part is combined 
with the heavy-quark part before it is dimcnsionally regularized at D = 4. Finally with 
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P 



OPE 



(*) = pP ert (a) + P ^ (s) + p^ 2 (a) + p^ Gs ^ (a) + P ^ (s) + p^ (a), 



p peit (s) = - — ^ J J -gj(l-a- (3)K(a,P)[r(mQ,8)-5mQrn s }r(m Q ,s) 4 , 



piSS){s) = S^^^S Si a ^{[( 2 -«-/3)m Q + (l-a-^)m s ]r(m Q , S ) 

/■"max ^ 

-3(a- a 2 + j3 - f3 2 )m s m 2 Q }r{m Ql s) 2 - m s / - [m Q - a(l - a)s} 3 }, 

p{ Ssf {s) = |Q^^ amOX da{ _ 2mQ[TO |_ a(1 „ a)s]+TOs[TO |_ 2a(1 _ a)s]}] 

p{g - s „.Gs ){s) = g |^ Q ^ r(mQiS){ _ 3mQ ( a + /3 -4a/3Mm Q , S ) 

+ m s a(3[l2m 2 Q - 7(a + /3)m Q - 5a/3s]} 

3m Q 
a(l — a) 

x [( m Q ~ 3m s )f(mQ, s) - 2m s mg(a + /?)], and 
f^OO = P"^ 1(1 - « - «%«{[(« 3 + ^K, + 4(« 4 + 



da[rriQ — a(l — a)]{— — ^[mg — a(l — a)s] + 2m s [5a(l — a)s — 9mg]}}, 



- 2m Q TO s (2a 2 + 3a£ + 2^ 2 )(3a 2 - 4a/3 + 3/3 2 )]r(m Q , s) - 4m s m Q (a + /3)(a 4 + /3 4 )}. 

It is defined as r(rriQ, s) = (a+f3)m,Q—a[3s and K(a, f3) = l+a—2a 2 +f3+2af3—2[3 2 . The integration limits 
are given by a min = (1 - ^1 - 4mg/s)/2, a maa; = (1 + ^1 - 4toq/s)/2, and /3 m i„ = amg/(sa - mg). 

III. NUMERICAL ANALYSIS AND DISCUSSIONS 

In this section, the sum rule ([7]) will be numerically simulated. The input parameters are taken as 
{qq) = -(0.23±0.03) 3 GeV 3 (ss) = 0.8 {qq), (gsa ■ Gs) = ml {as), to 2 , = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 , 



and {g 3 G 3 ) — 0.045 GeV 6 24, 28, 38|. For the quark masses, we employ the same values as Ref. [39j and 
references therein, which spanned by the running MS mass and the on-shell mass from QCD sum rule, 
with m c = 1.26 ~ 1.47 GeV, m b = 4.22 ~ 4.72 GeV, as well as m s = 114.5 ± 20.8 MeV. Complying with 
the standard criterion of sum rule analysis, the threshold so and Borel parameter M 2 are varied to find the 
stability window. It is well known that the fundamental assumption of the QCD sum rule is the principle 
of duality: it is assumed that there is an interval over which a hadron may be equivalently described at 
both the quark level and the hadron level. Therefore, the correlation function is evaluated in two different 
ways: at the quark level in terms of quark and gluon fields and at the hadronic level. If both sides of 
the sum rule were calculated to arbitrarily high accuracy, the matching of them would be independent of 
M 2 . Practically, however, both sides are represented imperfectly. On one hand, there are approximations 
in the OPE of the correlation functions and, on the other hand, there is a very complicated and largely 
unknown structure of the hadronic dispersion integrals in the phenomenological side. Thus, the extracted 
result is not completely independent of M 2 . The hope is that there exists a range of M 2 , in which the 
two sides have a good overlap and information on the resonance can be extracted. In practice, one can 
analyse the OPE convergence and the pole contribution to determine the allowed Borel window of M 2 : the 
lower limit constraint for M 2 is obtained by restricting that the perturbative contribution should be larger 
than the condensate contributions; the upper limit constraint is gained by the consideration that the pole 
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contribution should be larger than QCD continuum contribution. Meanwhile, the threshold parameter 
y/so characterizes the beginning of the continuum state. Thereby, it is not arbitrary but correlated to the 
energy of the next excited state with the same quantum number as the studied state. 

At first, we keep the values of the quark masses and condensates fixed at the central values. The 
comparison between pole and continuum contributions from sum rule ^ for Yj cs i for ^/sq = 5.2 GeV is 
shown in the left part of FIG. 1, and its OPE convergence by comparing the perturbative, quark condensate, 
four-quark condensate, mixed condensate, two-gluon condensate, and three-gluon condensate contributions 
is shown in the right one. Numerically, the ratio of perturbative contribution to the total OPE contribution 
at M 2 = 2.5 GeV 2 is nearly 60%, which is increasing with the M 2 to insure that perturbative contribution 
can dominate in the total OPE contribution when M 2 > 2.5 GeV 2 . On the other side, the relative pole 
contribution is approximate to 52% at M 2 = 3.2 GeV 2 and descending along with the M 2 to guarantee 
the pole contribution can dominate in the total contribution while M 2 < 3.2 GeV 2 . Thus, the region of 
M 2 for Y [cs] is taken as M 2 = 2.5 ~ 3.6 GeV 2 for y/s^ = 5.2 GeV. Similarly, the proper range of M 2 is 
gained as 2.5 - 3.0 GeV 2 for yfs^ = 5.0 GeV, and the range of M 2 is 2.5 - 3.6 GeV 2 for yfs^ = 5.4 GeV. 
We see also that for y/so = 4.9 GeV, the corresponding Borel parameter range is M 2 — 2.5 ~ 2.7 GeV 2 , 
which is very narrow as a working window. It is the main reason that ^/ sp < 4.9 GeV is not chosen 
here. In order to evaluate the uncertainty of results more conservatively [40j, we enlarge the variation of 
threshold parameter y/so for Y cs from 5.0 ~ 5.4 GeV to 5.0 ~ 5.7 GeV and we find the range of M 2 is 
2.5 ~ 3.8 GeV 2 for y/so — 5.7 GeV. In the chosen region, the mass result is not completely independent of 
M 2 since both sides of the sum rule are not calculated to arbitrarily high accuracy but have included some 
approximations, and that is just the reason by which the accuracy of QCD sum rule method is limited. 
Whereas, it is expected that the two sides have a good overlap and information on the resonance can be 
safely extracted in the chosen range of M 2 . The corresponding Borel curve to determine the mass of Y[ cs ] 
is exhibited in the left part of FIG. 3. We compute the average mass value of these working windows as 
4.69±0.29 GeV (the numerical error reflects the uncertainty due to variation of so and M 2 ). Up to now, we 
have kept the values of the quark masses and condensates at the central values. At last, we vary the quark 
masses as well as condensates and arrive at 4.69 ± 0.29 ± 0.07 GeV (the first error reflects the uncertainty 
due to variation of so and M 2 , and the second error resulted from the variation of QCD parameters) or 
4.69 ± 0.36 GeV in a concise form. 

For Yrf, s i , the comparison between pole and continuum contributions from sum rule (j6|) for y^so = 
11.8 GeV is shown in the left part of FIG. 2, and its OPE convergence by comparing different OPE 
contributions is shown in the right one. In detail, the perturbative contribution versus the total OPE 
contribution at M 2 = 7.5 GeV 2 is nearly 62%, and the relative pole contribution is approximate to 50% 
at M 2 = 9.0 GeV 2 . Thus, the region of M 2 is taken as M 2 = 7.5 - 9.0 GeV 2 for y/sE = 11.8 GeV. With 
the similar analysis, for y/so — H-6 GeV, the range is M 2 = 7.5 ~ 8.3 GeV 2 ; for y/so — 12.0 GeV, the 
range is M 2 = 7.5 ~ 9.5 GeV 2 . To evaluate the uncertainty of results more conservatively, we enlarge 
the variation of ^ horn 11.6 - 12.0 GeV to 11.6 - 12.3 GeV. For y/sE = 12.3 GeV, the range of M 2 
is 7.5 ~ 10.3 GeV 2 . The dependence on M 2 for the mass of Yf^i from sum rule (J7]) is shown in the right 
part of FIG. 3. For Y& s , We arrive at 11.19 ± 0.28 GeV (not including the variation of QCD parameters). 
Finally, we vary the quark masses as well as condensates and arrive at 11.19 ± 0.28 ± 0.21 GeV (the former 
error reflects the uncertainty due to variation of so and M 2 , and the latter error resulted from the variation 
of QCD parameters) or 11.19 ± 0.49 GeV in a concise form. 

With regard to the numerical results, some more discussions are given below. Numerically, the result 
4.69 ± 0.36 GeV for Yr cs i is in good agreement with the experimental value 4664 ± 11 ± 5 MeV for V(4660). 
However, its value is a bit higher than Y"(4260)'s mass even considering the uncertainty, which supports 
the P-wave [cs][cs] structure for V(4660) while disfavors the explanation of V(4260) as the P-wave [cs][cs] 
state. Note that some authors also assume that V(4260) could be a P-wave [cq] [cq] state In fact, 

we have calculated the mass of the P-wave [eg] [eg] to be 4.32 ± 0.20 GeV 4l|, which is compatible with 
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FIG. 1: In the left part, the solid line shows the relative pole contribution (the pole contribution divided by the total, 
pole plus continuum contribution) and the dashed line shows the relative continuum contribution from sum rule 
([6]) for yfso — 5.2 GeV for Yr cs i. The OPE convergence is shown by comparing the perturbative, quark condensate, 
four-quark condensate, mixed condensate, two-gluon condensate and three-gluon condensate contributions from 
sum rule ([6| for y/so = 5.2 GeV for Yj C3 ] in the right one. 




FIG. 2: In the left part, the solid line shows the relative pole contribution (the pole contribution divided by the total, 
pole plus continuum contribution) and the dashed line shows the relative continuum contribution from sum rule ((6]) 
for y/so = 11.8 GeV for Yfjgi. The OPE convergence is shown by comparing the perturbative, quark condensate, 
four-quark condensate, mixed condensate, two-gluon condensate and three-gluon condensate contributions from 
sum rule ([6]) for y/so = 11.8 GeV for Y[b s ] in the right one. 



the experimental data of Y(4360) and could support Y(4360)'s P-wave [cq][cq] structure. Barely from the 
value 4.32 ± 0.20 GeV, one could not completely exclude the possibility of Y(4260) as a P-wave [cq][cq] 
state since it is still in accord with the mass of Y(4260) in view of the uncertainty. Concerning the real 
nature of Y(4260), some further theoretical study and experimental verification are undoubtedly needed. 



IV. SUMMARY 

The QCD sum rule method has been employed to compute the mass of P-wave [cs][cs] tetraquark state 
Yr cs i, including contributions of operators up to dimension six in the OPE. The final result 4.69 ±0.36 GeV 
(4.69 ± 0.29 ± 0.07 GeV, where the first error reflects the uncertainty due to variation of sq and M 2 , 
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FIG. 3: In the left part, the dependence on M 2 for the mass of Yj cs j from sum rule (JT]) is shown. The continuum 
thresholds are taken as y/sE = 5.0 ~ 5.7 GeV. For y/sE = 5.0 GeV, the range of M 2 is 2.5 ~ 3.0 GeV 2 ; for 
y/sE = 5.2 GeV, the range of M 2 is 2.5 ~ 3.2 GeV 2 ; for y/sE = 5.4 GeV, the range of M 2 is 2.5 ~ 3.6 GeV 2 ; for 
y/so = 5.7 GeV, the range of M 2 is 2.5 ~ 3.8 GeV 2 . The dependence on M 2 for the mass of Yji, s i from sum rule (0 
is shown in the right one. The continuum thresholds are taken as y/so = 11.6 ~ 12.3 GeV. For = 11.6 GeV, 
the ranges of M 2 is 7.5 ~ 8.3 GeV 2 ; for y/IE = 11.8 GeV, the range of M 2 is 7.5 ~ 9.0 GeV 2 ; for yfs^ = 12.0 GeV, 
the range of M 2 is 7.5 ~ 9.5 GeV 2 ; for y/sE = 12.3 GeV, the range of M 2 is 7.5 ~ 10.3 GeV 2 . 



and the second error resulted from the variation of QCD parameters) for Y\ cs \ is well compatible with the 
experimental data of F(4660), which favors the P-wave tetraquark configuration for y(4660). Meanwhile, 
the result is higher than Y"(4260)'s mass, which is not consistent with assumption of Y"(4260) as the P-wave 
[cs] [cs] state. As a byproduct, the mass for the bottom counterpart Y[f, s ] has also been predicted, which is 
11.19 ±0.49 GeV (11. 19 ±0.28 ±0.21 GeV, where the former error reflects the uncertainty due to variation 
of so and M 2 , and the latter error resulted from the variation of QCD parameters) and expecting further 
experimental identification. 
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